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ABSTRACT 

We use the lag-luminosity relation to calculate self-consistently the redshifts, 
apparent peak bolometric luminosities L^, and isotropic energies i^iso for a large 
sample of BATSE gamma-ray bursts. We consider two different forms of the lag- 
luminosity relation; for both forms the median redshift for our burst database is 
1.6. We model the resulting Ei^o sample with power law and Gaussian probability 
distributions without redshift evolution, both of which are reasonable models. 
The power law model has an index of aE = 1-76 ± 0.05 (95% confidence), where 
p{Eiso) E^^'^ . The simple universal jet profile model suggested but did not 
require as = 2, and subsequent physically reasonable refinements to this model 
permit greater diversity in aE, as well as deviations from a power law; therefore 
our observed Eis^ probability distribution does not disprove the universal jet 
model. 

Subject headings: gamma-rays: bursts 

1. Introduction 



The major breakthroughs in the study of gamma-ray bursts of the past six years — most 
if not all bursts are cosmological, bursts do not have constant peak luminosity, the fireballs 
are beamed, many bursts are associated with supernovae — resulted from the intensive study 
of a relatively small number of bursts without regard for whether these bursts formed a well 
defined statistical sample. However, realizing that for any one burst we are only sampling 
an anisotropic radiation pattern from one direction, we now have to study distributions of 
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burst properties to reconstruct the bursts' appearance from all directions. Well-defined burst 
samples are required to derive these distributions. 

Unfortunately, we do not yet have a large sample of bursts with spectroscopic redshifts, 
which in most cases are required for calculating the intrinsic burst properties. The two 
dozen or so bursts with redshifts were detected by various detectors with different detection 
thresholds and energy sensitivities, and the follow-ups that determined the redshifts de- 
pended on the vagaries of weather and telescope availability. BATSE provided a large burst 
sample with well-understood thresholds, but without direct redshift determinations for most 
of these bursts. However, redshifts can be determined indirectly from the lag- luminosity 
(Norris, Marani & Bonnell 2000) or variability- luminosity (Fenimore & Ramirez-Ruiz 2000; 
Reichart et al. 2001) relations. Here we calculate self-consistently the redshifts for a large 
subset of the BATSE bursts using the lag-luminosity relation. To accommodate the prob- 
lematic burst GRB980425, which is significantly underluminous if the lag-luminosity relation 
is a single power law, Salmonson (2001) and Norris (2002) proposed that the relation should 
be a broken power law. This reduces the luminosity of many long-lag bursts, moving this 
population closer to the observer. 

We use the resulting redshifts to calculate both the peak bolometric luminosity Lb and 

the isotropic energy Eiso- Bursts are thought to radiate anisotropically, but we sample their 
radiation field in only one direction. Therefore Lb and i^iso are calculated from the observed 
fiux as if the emissions are isotropic; corrections for the anisotropy are based on models of 
the relativistic jets that emit the observed gamma rays (see Bloom, Frail & Kulkarni 2003). 
Lb and ii^iso are both bolometric quantities in the burst frame. is the maximum value of 
the luminosity while Eiso is the total energy emitted over the duration of the burst (without 
correcting for the anisotropic emission). We consider the isotropic energy Eiso to be more 
fundamental, and therefore model its probability distribution. 

Two models have been proposed for the structure of the jets. The uniform jet model 
(Frail et al. 2001; Bloom et al. 2003) assumes that all jets have a constant surface energy 
density e (energy per sohd angle across the jet) but differing opening angles 6o (the angle 
between the jet axis and the edge of the jet). The model makes no predictions for the 
energy probability distribution. On the other hand, in the universal jet profile model (Rossi, 
Lazzati, & Rees 2002; Zhang & Meszaros 2002) all jets have the same surface energy density 
e as a function of off-axis angle 6 (the angle from the jet axis), and thus the observed 
differences in Lb or Ei^o result from the angle 9^ between the jet axis and the line-of-sight. 
This model predicts that the energy probability distribution is a power law with index a^, 
where p(^iso) oc Er^"^. If e oc e^ then = 1 - 2/k. Rossi et al. (2002) and Zhang 
& Meszaros (2002) suggested k — —2, resulting in = 2, to reproduce the distributions 
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observed by Frail et al. (2001), although this was not a firm prediction. A Gaussian surface 
energy density profile results in a^; = 1 (Lloyd- Ronning, Dai & Zhang 2004). Thus a^; is 
of order 1-2. Lloyd-Ronning et al. (2004) showed that allowing the parameters of the jet 
profile to vary somewhat from burst to burst — the jet profile is now only quasi-universal — 
results both in variations in the value of a^; and deviations from a pure power law. While 
numerical modelling of hypernovae (Zhang, Woosley, & MacFadyen 2003) shows that the 
surface energy density of the outflows should vary with off-axis angle, it does not predict the 
profile; thus the profiles are fitted empirically to the data. Consequently, the observed E'lso 
probability distribution is not a powerful discriminant between the two jet structure models. 

This paper has two goals. First, we calculate the redshift of a large sample of bursts 
self-consistently from the lag-luminosity relation. With the redshift and the observed fiuxes 
and fiuences we calculate Lb and E'iso- Second, we model the Eiso distribution with two 
functional forms: power law and lognormal. We find the best fit values for the parameters 
of each functional form, and evaluate how well each functional form describes the data. We 
also compare the fit of the power law functional form to the predictions of the universal jet 
profile model. 

In these calculations we use a cosmology with Hq — 70 km s~^ Mpc~^, = 0.3 and 
Q,\ = 0.7. The notation p(a | b) means the probability of a given b. Lower case p denotes 
a probability density while upper case P (without a subscript) represents a cumulative 
probability. Note that a and b are propositions that arc true or false. A proposition can be 
a statement such as "the energy distribution can be described with a lognormal functional 
form" or a parameter value (equivalent to the statement "the parameter value is 5"). 

In §2 we present the methodology used in this study: calculating the burst rcdshifts 
(§2.1); fitting the energy probability distribution (§2.2); and using the cumulative probability 
to test the quality of the fit (§2.3). The implementation of this methodology is discussed in 
§3 and the results are provided by §4. Finally, our conclusions are in §5. 



2. Methodology 

2.1. Calculating Redshifts 

In the absence of a large number of measured spectroscopic redshifts, we use the lag- 
luminosity relation (Norris et al. 2000) to calculate peak luminosities Lb, and subsequently 
redshifts, for a large burst sample. Assume that for a set of bursts we have measured the lag 
To (s) between the light curves in two energy bands, the peak photon fiux Pph (ph s~^ cm~^) 
in the photon energy band El-Eu, and the spectral fit for the time period over which Pph is 
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measured. The spectral fit (A^ph(-E') — ph cm~^ keV~^) uses the "Band" functional form 
(Band et al. 1993) which is characterized by an asymptotic low energy power law with index 
a (A^ph(-E) oc -E"), a high energy power law with index /5 (iVph(£') oc and a characteristic 
energy between the two power laws; for /3 < —2 E^ is energy at the peak of the E'^N^n 
(energy per logarithmic frequency interval) light curve. This functional form is fitted to the 
observed count spectrum, which covers a limited energy range, but is extrapolated to high 
and low energy, as needed. The instrument team is assumed to have converted the observed 
peak count rate into the peak photon flux using the fitted spectral model and the response 
function. 

The lag tb in the burst's frame and the observed lag Tq, both between the same energy 
bands in their respective frames, are related by time dilation (which increases the lag) and 
spectral redshifting (which shifts the smaller lag at high energy into the observed band). 
The first effect is simple and universal — a factor of {1 -\- z)~^. The second depends on the 
burst's spectral evolution, and may vary from burst to burst. Thus 



We assume that g{z) is universal, and as a working assumption we use g{z) = + . Time 
dilation contributes —1 to c^-, while the redshifting of temporal structure with a smaller lag 
from higher energy contributes a positive constant (e.g., ~ 1/3; Fenimore et al. 1995 found 
~ 0.4). 

Empirically the lag has been related to the apparent bolometric peak luminosity Lb 
(erg s~^; Norris et al. 2000) 



The observed peak bolometric energy flux Fb (erg cm ^ s ^) is related to the bolometric 
peak luminosity Lb through the redshift, and an assumed cosmology: 



where is the luminosity distance. Note that Lb is the "isotropic" peak bolometric 
luminosity, the peak luminosity if the observed flux were beamed in all directions. If the flux 
is actually beamed into a solid angle AQ, then the actual peak luminosity is only Afl/An of 
Lb- We deflne 



Tb = rog{z) 



(1) 



Lb = Q(ts) 



(2) 



Fb = Lb/IAttDI] 



(3) 



(E) 



(4) 



where 



and 




(5) 



(6) 
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A^ph(-E') is the photon spectrum. The photon spectrum energy flux Fb results from a fit to 
the observed spectrum which is available over a limited energy band from a time bin that 
includes the peak of the light curve. Thus the peak bolometric energy flux Fb is calculated 
by extrapolating Npi^{E) in eq. 6 to high and low energies. 

The resulting implicit equation for redshift, 

p _ Q {rog{z)) 

must be solved for each burst. The inputs are tq, Pph and (E) (which is calculated from 
the spectral fit). The functional form of the lag- luminosity relation Q is calibrated from the 
small set of bursts for which z is known. We have an assumed functional form for g{z) with 
one unknown constant c^. This constant can be calculated using the dependence of r on 
energy bands (Fenimore et al. 1995). 

These equations can be evaluated for the sensitivity of calculated quantities such as 
z and Lb on observables such as Pph and (E) (which, in turn, depends on the spectral 
parameters). Eq. 7 shows that z — (f){Ppi^{E) , tq) while eqs. 3-4 can be combined to give 
Lb — Pph{E)ip{z), where and ip are functions defined only for this sensitivity analysis and 
need not be derived explicitly. The quantities Pph and (E) appear in these expressions only 
in the product Pph(£'), and consequently errors in (E), resulting from uncertainties in the 
spectral parameters, are equivalent to errors in Pp^. As shown by figure 1, these equations 
give curves in z-Lb space parameterized by Tq; larger values of the product Pph(£') lie at 
larger values of z and on the appropriate curve. 

If the spectrum is fairly fiat, e.g., a = —1 and f3 = —2, then {E) does not vary by a large 
factor as Ep varies. The segments plotted on figure 1 arc for Ep ranging between 5 kcV and 
2 McV for « = —1 and f3 = —2 while holding Pph fixed. However, if the spectrum is more 
peaked, i.e., is harder at low energies (larger a) or softer at high energy (smaller f3), then 
the (E) range is larger. This is shown on figure 1 by the diamonds, triangles and asterisks 
that mark the range that results as Ep varies for different sets of a and p. As can be seen, 
cavalier treatment of the spectrum of the peak flux can result in large errors in z and Lb- 

Solving eq. 7 provides z for each burst (note that Lb does not appear explicitly in eq. 
7). With Pph and the energy fluence (erg cm^^) we then calculate Lb and E'iso- The result is 
a large database of observables (lags, Pph, energy fiuences, spectral parameters) and derived 
quantities {z, Lb, and Eiso) with which we can address different questions about the burst 
population. Here we study the distribution of intensity measures such as Lb and Eiso- 
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2.2. The Intensity Probability Distribution 



We use a Bayesian analysis to investigate intensity probability distributions, although 
we argue that a frequentist analysis will give the same result; Band (2001) provides a more 
complete exposition of the methodology. In this study we consider the isotropic energy i^iso 
as the relevant intensity measure; the same formulae apply substituting Lb for E^so- 

Our burst database provides the isotropic energy -Eiso, the peak photon flux Pph and the 
redshift z for each burst. The threshold peak flux Pph,min is known, from which the threshold 
isotropic energy -Eiso,min for each burst is calculated: £'iso,min = -E-iso -Pph,min/-Pph- The energy 
probability distribution function is the normalized probability distribution p(-E'iso | CLji Mj, I) 
where d} is the set of parameters that characterize the jth model distribution function 
represented by Mj, and / specifies general assumptions about our calculation. Thus Mj 
states that we are using the jth functional form, which has parameters dj. There are always 
additional assumptions upon which the calculation rests, such as the underlying cosmology 
or the detector calibration. The validity of our analysis depends on these assumptions, which 
we represent with the proposition /. If the distribution evolves with redshift, then z should 
be included in the list of 'givens' in the probability distribution; here we do not include a 
redshift dependence. Thus we do not consider redshift evolution in our analysis. Because we 
work with probability distributions that are normalized at each redshift, our results remain 
valid for evolution in the burst rate per comoving volume. 

In our calculations we use power law and lognormal energy probability distributions. In 
the power law case the proposition Mj = Mpi is the statement that the functional form is a 
power law, and the parameters Uj are the power law index and the low energy cutoff E2: 



P{Eiso I OiE, E2, Mpi, /) 



oc Er^- for 



Eiso > E2 



(8) 



For a finite number of bursts the distribution must have a low energy cutoff £'2 if > 1. As 
we will see, the fit to the power law model is relatively insensitive to the low energy cutoff 
E2. Therefore our results are vahd for redshift evolution in the energy scale (i.e., in E2), but 
not in the power law index. 

In the lognormal case Mj — Min states that the functional form is a lognormal, and the 
parameters are the energy centroid E'iso.cen and the logarithmic width (Tg: 



p{Ei 



, E- 

ISO -"^isoxen 



, (7e, Min, /) = 



1 



27r(7E 



exp 



(In (Ei^) - ln(£;iso,cen))^ 



2al 



(9) 



For the ith burst we have £^iso,i with threshold £^iso,min,i, the set of which constitutes the 
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data D. The observed £'iso,i are not drawn from p(£'iso | Mj, I) but from 

where H{x) is the Heaviside function (1 for positive x and for negative x); E^go is drawn 
from the observable part of the energy distribution function. 

The probabihty of obtaining the data D given the model is the 'likehhood' 

Nb 

Aj = p{D I Uj, Mj, I) = Y[p{Ei,o,i I £^iso,min,i, %, Mj, I) (11) 

i=l 

where the product is over the Nb bursts in our database. In the "frequentist" framework 
best-fit parameters are typically found by maximizing through varying the parameters. 
Plotting Aj as a function of (Tj around the maximum reveals the range of acceptable dj 
values. 

However the Bayesian analysis is based on p{dj \ D, Mj, I), the posterior probability for 
the parameters, that is, the probability that the particular set of dj values is correct given 
the data. Bayes theorem gives 



P(D\ 




.Mj.I)p{aj\ 




jdajp{D\ 


\d'j,Mj,I)p{d'j 


Mj, I) 



The factor p(d} | Mj, I) is the prior for dj, constraints on the parameters based on information 
available before the new data were acquired. The factor in the denominator is a normalizing 
constant. Thus for parameter determination the Bayesian approach explicitly factors in 
additional information and constraints (e.g., energies must be non-negative and finite). The 
expectation value of the parameters is 

(a-)=/<ia-«-p(«-|AM„/) . (13) 

If Sj = p{D I d'j, Mj, I)p{dj I Mj, I), the numerator in eq. 12, is sharply peaked, then the ex- 
pectation value of the parameters occurs at the peak of Ej. However the posterior probability 
p{dj I D, Mj, I) (or Ej) is required to determine the acceptable parameter range. 

As was discussed in Band (2001), reasonable priors for the parameters of both the power 
law and lognormal distributions are constants (when the logarithm of an energy is treated 
as the parameter rather than the energy), and therefore the frequentist likelihood Aj and 
the Bayesian posterior p{dj \ D, Mj, I) are proportional to each other. Consequently, while 
we favor the Bayesian approach, here the frequentist and Bayesian analyses are the same. 
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We calculate as a function of the parameters d}. The maximum of Aj gives the best fit 
parameter values while the width of shows the range of acceptable values, and correlations 
between the parameter values. 

Note that the methodology we use here only provides the normalized probability dis- 
tribution, not the normalization of the distribution (i.e., the burst rate per volume). The 
normalization requires the history of £'iso,min during the mission (i.e., not only for the detected 
bursts) . 



2.3. The Cumulative Probability 

The Bayesian approach does not provide a goodness-of-fit statistic. However a frequen- 
tist statistic can be derived. For each burst the cumulative probability is 

POD 

-^(-^ -^iso,i I -^iso,min,i) 

d}, Mj, I) = I p{Eiso I -Eiso,min,i, (^j, Mj, I) dEigo . (14) 

If the assumed energy distribution function is an acceptable characterization of the observa- 
tions (which would be the case if the model Mj is correct) and all the assumptions (included 
in the proposition /) are valid (e.g., the cosmological model is correct), then the cumulative 
probabilities P{> £'iso,i) for each burst should be uniformly distributed between and 1, 
and have an average value of (P(> -Eiso.i)) = 1/2 ± {12Nb)~^^'^ for the Nb bursts in the 
sample. Note that {12Nb)^^^'^ is the expected statistical variance resulting from the size of 
the sample, and does not take into account systematic errors. 



3. Implementation 
3.1. Datasets 

We start with a database for 1438 BATSE bursts that includes the lags and their 
uncertainties, the peak flux Pph over the 50-300 keV band on the 256 ms timescale, the 
burst duration Tgo, and hardness ratios among the 4 BATSE energy channels (30-50, 50- 
100, 100-300, and 300-2000 keV). Of these 1438 bursts, 1218 have positive lags. 

To calculate the average energy (E) (see eq. 4) we use the parameters of the "Band" 
spectrum fits by Mallozzi et al. (1998) to the 16 channel BATSE "CONT" count spectra 
accumulated over the peak flux time interval (usually 2.048 s) for the 580 bursts in our 
database that are also in the Mallozzi et al. database. For the 858 bursts that are not in 
the Mallozzi et al. database we assume the average low and high energy spectral indices of 
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a — —0.8 and (5 — —2.3 found by Preece et al. (2000). Because Ep^ the energy of the peak 
of the E'^N^^{E) (oc vf^) curve, and HR32, the 100-300 keV to 50-100 keV hardness ratio, 
are strongly correlated, we use the empirical relation Ep = 240 HR32 keV (for HR32 < 2.25) 
for the bursts without spectral fits. In §2.1 we discussed the sensitivity of z to errors in the 
spectrum; the gaps in the available spectral fits introduces a systematic uncertainty into our 
analysis. Figure 2 shows the resulting scatter plot for Ep vs. HR32. 

3.2. The Lag-Luminosity Relation 

Based on less than a dozen bursts, Norris et al. (2000) found that Lb oc t^^'^^. However, 
this single component relationship predicts a much greater Lb than observed for GRB980425, 
the burst that appears to coincide with the supernova SN1998bw. By introducing a break 
in the lag- luminosity relation Salmonson (2001) and Norris (2002) were able to include 
GRB980425. The resulting lag-luminosity relation is: 

L51 = 2.18(rB/0.35 s)'=^, cl = -1.15 for tb < 0.35 s, - 4.7 for tb > 0.35 s. (15) 

where L51 = L^/IO^^ erg s~^. In our calculations we use both the single (i.e., cl — —1.15 
for all Tb) and two component lag-luminosity relations. 

4. Results 

4.1. Lag-Luminosity Relation and Redshift 

We first investigate the lag- luminosity relation and then use it to construct a database 
of burst redshifts. Lb and E'iso- Salmonson (2001) and Norris (2002) introduced a break 
in the simple power law relation to include the long lag, very low luminosity GRB980425. 
As a result of this break other long lag bursts in the BATSE database are assigned low 
luminosities, and consequently small distances; Norris (2002) suggested that there may be 
a population of nearby low luminosity bursts. Figure 3 demonstrates this: the solid line 
is the cumulative distribution of burst redshifts assuming a single power law lag-luminosity 
relation, while the dashed line shows the distribution for the broken power law relation. 
Introducing a break in the power law shifts the low redshift bursts closer (i.e., to lower 
redshifts). Note this break was introduced only to include GRB980425 in the lag- luminosity 
relation; if GRB980425 was not associated with SN1998bw or if GRB980425 was anomalous, 
then this break is unnecessary and unsupported by other data. 

Figure 3 also shows that there are few bursts with 2; > 10. These bursts may indeed 
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be at high redshift; finding such high redshift bursts is one of the goals of the Swift mission. 
However, on theoretical grounds bursts with z > 17 are not expected. Given the dispersion 
in the lag-luminosity relation and the errors in determining tq, it is not surprising that some 
bursts are assigned redshifts as high as z = 65; large redshifts result from high luminosities 
for bursts with very small lags, which are particularly difficult to measure accurately. Thus 
the lag-luminosity relation does not give unphysical results at the high luminosity end. Norris 
(2002) placed an upper limit on the luminosity when he found a large number of very high 
redshift bursts. 

We use the simple single power law lag-luminosity relation without any cutoffs or limits 
in the following analysis. 

For our database the median redshift is = 1.58. As can be seen from Figure 3, 
this redshift is greater than the redshift range where the form of the lag-luminosity relation 
makes a difference. 

Figure 4 shows the distribution of we calculate (the scatter plot) along with the 
threshold of at different redshifts (the lines). In calculating the threshold values of Lb 
we assume a constant Pph.min = 0.3 ph cm~^ s~^ and vary (E) (see eq. 4). Factor of ~ 3 
differences in (E) result in significant over- or under-predictions of the threshold. 

Thus we have accomplished the first goal of this paper. The resulting database is 
provided online at http://cossc.gsfc.nasa.gov/analysis/lags/ . 

4.2. The ii^iso Probability Distribution 

We focus on the probability distribution functions of the bolometric energy. To perform 
quantitative estimates of distributions as a function of redshift (e.g., of the burst rate) the 
threshold for the database must be known, and therefore studies have often focused on the 
peak luminosity Lb', Lb is more closely related to the peak fiux Pph, which has a relatively 
sharp instrumental threshold Pph,min- However, to study probability distributions of intrinsic 
burst quantities the threshold for that quantity for each observed burst will suffice. Thus 
our database is sufficient to study both the isotropic energy E^go and the peak luminosity 
Lb probability distributions. 

Figure 5 shows the scatter plot of the isotropic energy Eiso- The region just above 
the threshold is underpopulated, suggesting that the thresholds are underestimated for this 
database. Therefore, we raise the threshold to Pph,min = 0.5 ph cm~^ s~^, raising the median 
redshift from 1.58 to 1.62, and decreasing the number of bursts to 1054. 
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4-2.1. Power Law Eiso Distribution 

This functional form is relevant because it is predicted by the universal jet profile model. 
Figure 6 shows the resulting likelihood surface, which peaks at = 1.76 and log(£'2) = 50.2. 
This value of E2 is the smallest value of ii^iso in our sample; burst samples that extend to 
fainter bursts will probably have smaller observed values of E'iso, and consequently E2 is 
most likely smaller. The contour plot indicates that the determination of ccg is independent 
of the value of E2. Since the likehhood peaks at the same q;^; at any given £'2, our result 
holds even if the energy scale evolves with redshift. The 95% confidence region centered on 
ue — 1.76 has a half width of 0.05. 

As discussed in §1, the universal jet profile model, especially with physically reasonable 
refinements, does not predict definitively a value of ckb that can then be used to falsify this 
model. Originally as = 2 was suggested (Rossi et al. 2002, Zhang & Meszaros 2002) for a 
surface energy density that is a power law in the jet off- axis angle, while a Gaussian surface 
energy density would have «£■ = 1 (Lloyd- Ronning et al. 2004). Consequently our value of 
aE = 1.76±0.05 falls in the range of expected aE values. Thus our work does not distinguish 
between the universal jet profile model and the uniform jet model (which does not predict 
an energy probability distribution). 

Since the value of is relevant to jet models, consideration of the possible systematic 
errors is warranted. Our method for estimating the energy probability distribution considers 
where each measured £"180 falls between £iso,min for that burst and infinity. If ^iso.min is 
underestimated, then the low ^igo portion of the probability distribution will be under- 
represented, causing the distribution to be shifted to higher energy, in this case resulting in 
a smaller value of ag. Indeed, when we use Pph,min=0.1, 0.3, 0.5, 1 and 2 ph cm~^ s~^ we find 
q;e=1.125, 1.6, 1.76, 1.9 and 2. Next, uncertainties in the measured value of E^so that are 
symmetric to higher and lower values will cause a net 'diffusion' towards higher i^iso because 
there are more bursts below than above any given value of E^so- This effect also decreases 
aE- Note that small lags are more difficult to measure, resulting in greater uncertainties for 
large Lb; since E^so is correlated with Lb (but is not strictly proportional), on average the 
uncertainty increases for larger Eiso- Therefore, the true value of is most likely greater 
than the value we measured. 

The likelihood alone does not indicate whether the power law model is a good description 
of energy distribution. The average of the cumulative probability (§2.3) should be {P{> 
Eiso)) = 1/2 ± (12iVB)-V2 for bursts. Figure 7 shows that the cumulative P{> Eiso) for 
-Pph,min = 0.5 ph cm~^ s~^ is very close to the straight line that is expected if the distribution 
function is a good description of the data. We calculate {P{> -Eiso)) — 0.4642 with an 
expected statistical variance of 0.0089 {Nb — 1054). Thus {P{> -Eiso)) differs from the 
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expected value of 1/2 by 4(T using the statistical variance only. However, since the likely 
systematic uncertainties (e.g., in the value of i-*ph,mm, in the determination of the lags, in 
the corrections for the redshifting of high energy light curves) are likely to be considerable, 
this value of {P{> Eiso)), and the proximity of the P{> -Eiso) distribution to a straight line 
(Figure 7), indicates that a power law is a good representation of the distribution. 



4-2.2. Lognormal E^so Distribution 

The second model distribution we consider is a lognormal energy probability distribu- 
tion; the parameters are the energy centroid -Eiso.cen and the logarithmic width aE of the 
distribution. For this distribution the likelihood (relevant to a frequentist analysis) and the 
posterior distribution as a function of lni?iso,ccn and (relevant to a Baycsian analysis) arc 
the same if the prior is constant in lni?iso,ccn and (see Band 2001). Figure 8 shows the 
likehhood surface; the peak occurs at -Eiso,cen = 2.8 x lO^'^ ergs and erg = 2.7. However, these 
parameters are highly correlated since a broader distribution (larger ue) can compensate for 
a smaller central energy -Bisccen- The figure shows that 95% of the probabihty distribution 
is in the range -Biso.cen =(0.03-1) x 10^^ erg and ge = 2.3-3.2. As can be seen from figure 5, 
-E'iso.cen = 2.8 X 10^° crgs is at the lower end of the distribution of measured i?iso,min! and thus 
the data are insufficient to determine whether the probability density does indeed decrease 

below -Eiso.cen- 

For Pph,mm = 0.5 ph cm-2 (P(> Ei^o)) = 0.4821 ± 0.0089 (A^^ = 1054), which is 
consistent with (P(> -Eiso)) = 1/2 at the 2a level. Note that for Pph,min — 0.3 ph cm~^ s~^ 
(P(> Piso)) = 0.4598 ± 0.0085 {Nb = 1162), which differs from the expected value by 
nearly 5o" using the statistical variance only. As we argued above, there are undoubtedly 
significant systematic uncertainties in addition to the statistical variance. Figure 9 shows 
the distribution of the cumulative probability for Pph,min = 0.5 ph cm~^ s~^; as can seen, 
the observed cumulative probability distribution is very close to the expected distribution. 
Again, considering the systematic uncertainties in our calculation, the lognormal distribution 
is consistent with the data. 



4.2.3. Comparison to Previous Probability Distribution Calculations 

In this study we apply the same methodology presented in Band (2001) to a new burst 
data set. In both studies the E'iso distributions are modelled with lognormal and power law 
functional forms, although Band (2001) applied a high energy cutoff to the power law function 
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because in some cases < 1. Band (2001) used three data sets: a) 'B9' — 9 bursts with 
spectroscopic redshifts (i.e., the rcdshifts were measured from emission or absorption hues) 
and fitted BATSE spectra; b) 'C17' — 17 bursts with spectroscopic redshifts and spectral 
information from a variety of sources; and c) 'F220' — 220 bursts with redshifts derived 
from the variabihty-luminosity relation (Fenimore & Ramirez- Ruiz 2000). Table 1 presents 
the resulting parameters; 90% uncertainties are indicated. The variables of the lognormal 
distribution are highly correlated, resulting is large uncertainty ranges for both parameters. 
We find that as the database size increases, it samples smaller Eiso, as can be seen from the 
values of E2. Perhaps as a consequence, £'iso,cen decreases and as increases as the database 
size increases. 

The distributions for these different samples are discrepant, which may result from 
systematic difficulties with the burst samples or incorrect assumed probability distribution 
functional forms. The variabihty-luminosity relation used for the F220 sample and the lag- 
luminosity relation used in this paper were both calibrated with only a few bursts, and the 
validity of these relations must be confirmed by a larger burst sample. The B9 and CI 7 
samples use spectroscopic redshifts, and thus are affected by selection effects in detecting 
such redshifts. Band (2001) used the Pph,min for the BATSE detection, but the true Pph,mm for 
measuring the redshift was undoubtedly significantly greater: the bursts in this sample were 
usually detected and rapidly locahzed by a less sensitive detector (e.g., Beppo-SAX); and the 
burst's afterglow had to be sufficiently bright for an improved localization that warranted 
follow-up spectroscopic observations. Thus the B9 and C17 samples were seriously flawed. 
Note that the discrepancy between the F220 sample and ours is not very great, particularly 
for the power law distribution. 

The systematic trends in the fitted parameters of the probability distribution function 
with the burst sample's £^iso range suggest that a simple power law or a lognormal form are 
not the correct functional form. For example, Schaefer, Deng, & Band (2001) and Lloyd- 
Ronning, Fryer & Ramirez- Ruiz (2002) find that the luminosity function for is a broken 
power law. 

Comparison to other studies is more difficult, particularly since we study the probability 
distribution of E[so while others (Fenimore & Ramirez- Ruiz 2000; Schaefer ct al. 2001; 
Schmidt 2001; Lloyd- Ronning et al. 2002; Norris 2002) studied the luminosity function of 
Lb- In addition, these studies use different definitions of the luminosity, e.g., for BATSE's 50- 
300 keV trigger band (Fenimore & Ramirez-Ruiz 2000; Lloyd-Ronning 2002) or bolometric 
[E = 0.025-20 MeV for Norris 2002), with transformations using a variety of spectra (e.g., 
an E~^ power law for Schaefer et al. 2001 and Schmidt 2001, a broken power law for Norris 
2002 and a 'Band' spectrum with fixed parameters for all bursts for Fenimore & Ramirez- 
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Ruiz 2000). Consequently, a more detailed comparison between distributions is beyond the 
scope of our study. 

5. Summciry 

In this paper we have two objectives. First, we use the lag-luminosity relation to cal- 
culate self-consistently the redshifts for 1218 BATSE bursts. For the bursts without the 
spectral parameters required by the calculation we use average low and high energy spectral 
indices, and a peak energy Ep derived from the hardness ratio. We find that the redshift 
is quite sensitive to the spectral parameters, particularly if the spectrum is sharply peaked. 
We use both single power law and broken power law lag-luminosity relations, and find that 
the broken power law relation does indeed predict a population of low luminosity, nearby 
bursts. For both forms of the relation the median redshift is 1.58. 

We use the redshifts to calculate the apparent peak bolometric flux and the isotropic 
energy, both assuming that the bursts radiate isotropically. 

Second, we fit two functional forms to the distribution of the isotropic energy. We find 
that our burst data can be fit by a power law energy distribution with aE = 1.76 ± 0.05 
(95% confidence); considering the likely systematic uncertainties in addition to the statistical 
variance, the power law distribution is probably a good description of the data. This value of 
a^; is in the acceptable range for the universal jet profile model, and therefore, our work does 
not distinguish between the current jet structure models. A lognormal energy distribution 
also describes the data; the data permit a smaller average energy if the distribution is wider. 
More faint bursts will bound the lower end of this distribution. 

We thank the referee for helpful comments on improving the text of our paper. We also 
thank V. Avila-Reese for pointing out an error in the text. 
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Fig. 1. — Variation of bolometric peak luminosity as a function of redshift as Ep varies. 
The sohd curves, labelled by the observed lag tq, are for a = —1 and P — —2 while Ep 
varies between 5 keV and 2 MeV. Also shown on the tq = 0.1 curve and its extrapolation are 
the maximum and minimum redshifts for: triangle — a — —1, f3 — —2; diamond — a — — 1, 
P — —3; and asterisk — a — —1/2, P — —2. In all cases Pph = 1 ph cm" 
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Fig. 2.— Hardness ratio HR32 (100-300 keV vs. 50-100 keV) as a function of Ep. Two 
populations are evident. First, spectral fits were not available for the bursts that fall on an 
empirical Ep oc HR32 relation. Second, bursts with spectral fits are dispersed around this 
Ep oc HR32 relation. 
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Fig. 3. — Cumulative distribution of bursts by redshift. Solid line — one component lag- 
luminosity relation. Dashed line — two component relation. Both distributions have the 
same median redshift of 1.6. 
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Redshift 

Fig. 4. — Bolometric peak luminosity Lb vs. redshift z for the single component lag- 
luminosity relation. The lines represent the threshold peak luminosity for a threshold peak 
flux of Pph,min = 0.3 ph cm-2 s-^ in the 50-300 keV band, a = -1, /3 = -2, and {E) = 30 (2 
dots-dashed hne), 100 (dot-dashed), 300 (solid) and 1000 keV (dashed), where {E) is defined 
by eq. 4. 



-20- 




Fig. 5. — Scatter plot of E'iso vs. the -Biscmin- On the dotted hne Eiso — -Eiso,min- The 
paucity of bursts just above the dotted hne suggests that Pph,min should be raised from 0.3 
to 0.5 ph cm~^ s~^. 
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Fig. 6. — Contour plot of the likelihood function for a power law energy probability distri- 
bution. A threshold peak flux of Pph,min = 0.5 ph cm^^ s"^ was imposed. The power law 
has an index of aE and a low energy cutoff of E2. The contour levels contain 0.68, 0.9, 0.95, 
0.99 and 0.999 of the integrated probability. 
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Fig. 7. — Cumulative distribution of the probability of the observed isotropic energies for a 
power law energy probability distribution. A threshold peak flux of Pph,min = 0.5 ph cm~^ s~^ 
was imposed. 
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Fig. 8. — Contour plot of the likelihood function for a lognormal energy probability distri- 
bution. A threshold peak flux of Pph.min = 0.5 ph cm~^ s~^ was imposed. The distribution 
is characterized by a centroid energy and a logarithmic width. The contour levels contain 
0.68, 0.9, 0.95, 0.99 and 0.999 of the integrated probability. 
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Fig. 9. — Cumulative distribution of the probability of the observed isotropic energies for a 
lognormal energy probability distribution. A threshold peak flux of Pph,min = 0.5 ph cm~^ 
was imposed. 
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Table 1. Probability Distribution Parameters for Different Burst Samples 



Quantity 


B9^ 


C17'' 


F220^ 


LL1054^ 


-•^isOjCen 


130^}^° 

n 74+0.46 
'-'■'^-0.34 


CO +48 

0.96l°J? 


12+^^ 

-■-^-10 

-■-•^-0.2 
1 01+0.13 
-'-O-'--0.11 


9 7+0.4 
^- ' -0.3 

"J-0.07 


eJ^ 


1.6 


0.55 


0.12 


0.16 




1440 


1460 


5000 





^Sample of 9 BATSE bursts with spectroscopic redshifts and fitted spectra; analyzed in 
Band (2001). The threshold -E'iso,min was significantly underestimated. 

•^Sample of 17 bursts with spectroscopic redshifts; analyzed in Band (2001). The threshold 
-E'iso,mm was Significantly underestimated. 

■^Sample of 220 bursts with redshifts derived from variability redshifts; analyzed in Band 
(2001). 

•^The sample used in this work, with Pph,min = 0.5 ph cm~^ s~^. 

^The central energy of the lognormal distribution, in units of 10^^ erg. 

^The logarithmic width (in units of the energy's natural logarithm) for the lognormal 
distribution. 

^Thc index of the power law distribution, p{Eiso) oc E^^'^ . 

^The low energy cutoff of the power law distribution, in units of 10^^ erg. This energy is 
the lowest £^iso,min for the sample. 

'The high energy cutoff of the power law distribution, in units of 10^^ erg. 



